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1. Introduction 

A knowledge of the "scatter" ^or diffusion) of Individual hydrometeors 
(rain, snow, hall) is Indlspenst ble in the study of the effects caused by a 
population of hydrometeors on the propagation of radiowaves at frequencies 
above 10 GHz. 

I 

As a matter of fact all of the phenomena which influence propagation at 
such frequencies (attenuation, depolarization, scintillation, etc.) can be 
intensified by parameters directly derived from a solution of Individual 
/ 3 catt(?.r, naturally in addition to the meteorological elements which character- 
ize the physical medium. 

This paper will go into detail on the diffusion caused by rainy precipi- 
tation; the physical models which will be taken into consideration are rain 
composed of spherical drops, rain composed of drops in an ellipsoidal form 
with axes of rotational symmetry arranged along the vertical line of a 
generic reference point. 

* 

Since the most general type of electromagnetic wave polarization, includ- 
ing circular polarization and rectangular polarization as special cases, 
is elliptical polarization, there is a certain amount of interest in a 
detailed study, as complete as possible, of the effects of diffusion on tropo- 
spheric propagation of electromagnetic waves in the presence of rain alone 
with generic polarization, specified for the cases of interest to us within 
the "range" of frequency from 10 to 100 GHz. 

The case of earth-space propagation will not be treated in detail in the 
present work, although for the sake of completeness there will be a report of 
attenuation and phase displacement results obtained for some possible linkg 
with axes of propagation Inclined with respect to the vertical line in the 
case of rectilinear polarization alone. 
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2. Rain Composed of Spherical Drops: G. Mie's Algorithm 

The diffusion due to drops of a spherical shape will be treated briefly 
in this paragraph for the sole purpose of greater clarity and completeness of 
the results reported in this paper. 

G. Mle has developed an algorithm [B.l] to evaluate the individual scatter 
of spheres of any material whatsoever at any frequency whatsoever. 

For our specific interests the theory has been developed for spherical 
drops of rain at the F.U.B, [Fondazione Ugo Bordoni] [B.2]; the scatter 
functions in the "forward" direction have been evaluated from; 


oo 

C ' .^aln + Dn'\ (1) 

.2 ' ^ ' 

and the corresponding sections of extinction and phase displacement [B.2]: 

2 ^ 

Rt ^;'2(,2n*l)(a„ + b„) ( 2 ) 

In (2) D is the diameter (in ram) of the drops of rain taken into consideration, 
\ is the wavelength (in cm) , and a^ and b^ are the so-called coefficients 
of Mie. 

The specific attenuation and phase displacement have been determined by 
using the extinction and phase section [B.2]; 


As(d^m) = 4.34-10* (Ot)- N(p:) 

F 5 (y*di/Km) = ^ ^ C* (Di) • NCr>0 


(3) 


3 


i 


Every term In the preceding summations is the product of the extinction 

(or phase) section relative to drops of diameter by the number of drops 
3 

contained in 1 m with diameters in the field of D with AD«0.5 mm. 
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Tlie number of spherical drops N(D^) is found by means of the equation; 


• X w J 4 ' 


P • 


C* * I 


(4) 


in which R(nun/h) is the intensity of precipitation, M(D^) [B.3J Is the fraction 
of the total volume of rain due to drops with diamet-^rs in the field of 

I relative to the predetermined value of R, V(Dl) in m/sec is the 

X 2 

terminal falling velocity of the drops with diameter Di [B.4], The N(D^) 

presented in Table 2 were calculated from the values reported in Table 1. 

Assuming Ray's water refraction index [B.5] at a drop temperature of 20“G, 

the values obtained from A (dB/Km) and E (®/Km) for 19 frequencies between 

S Q 

10 and 100 GHz are presented in Tables 3-11 and 12-20 as a function of R. 


3. Rain Composed of Ellipsoidal Drops: T. Oguchi's Algorithm 

In the case of simultaneous transmission of two electromagnetic waves 
at the same frequency with one polarized rectilinear lead in the vertical 
direction and the other in the horizontal direction, experimental findings 
have demonstrated the existence of a difference in attenuation and in phase 
displacement between the two polarization directions [B.6]. 

The phenontenon cannot be explained with the simple physical model of 
spherical drops. 

For this reason it has been assumed that the real form of the raindrops, 
rather than being spherical, can be assimilated to an ellipsoid of rotational 
symmetry with axes arranged along the vertical line. 
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31iis then entails the need of going Into deeper studv of the scatter of 
an individual ellipsoid. Such a study has been carried out by T. Oguchi [B.7] 
by way of three different algorithms. Tlie one defined as the "Point Matching 
Technique" has been implemented at F.U.B., since it offers the advantage, over 
the other two, of being applicable to electromagnetic diffusors of any form 
whatsoever provided that they have rotational symmetry. 

This nusthod determines the coefficients of a series development of spherical 
harmonics analogous in principle to those of Mie, with the difference that, 
since in this case the diffusor is not spherical, the variable spherical 
harmonics are not decoupled: i.e., a definite spherical harmonic present in 
the incident wave is not limited to stimulating only the corresponding spherical 
harmonic in the diffuse wave, but stimulates a v^hole set of them coupled with 
it, in addition to stimulating spherical harmonics in the diffuse wave without 
any corresponding ones existing in the incident wave. 

Actually, the generic component of the electromagnetic field can be 
expressed as a function of the spherical harmonics which satisfy Maxwell's- 
equation [B.8j; all of the harmonics can be regrouped into two families, a and 
b, each with two degrees of freedom, m and n. 

The "point matching technique" is based on the finding that the special 
rotational symmetry of the problem guarantees that all of the harmonics of a 
different index m are decoupled from one another. This means that the harmonics 
can be regrouped into families characterized by the same index m, including 
harmonics coupled to one another by the various values of n. 

All of this considerably simplifies the problem, and in fact if N is the 
number of harmonics to be given consideration, it turns out that the two 
Indices, m and n, vary from one another: 
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a « SB < W 


(• ^o) 


■%<>,< 


Each m family therefore has N-wH unknown harmonic amplitudes, and therefore 

i 

requires N-m+1 independent equations [B.8]. 

Since there are four successive equations for the contour of the drop 
(two tangential components of E and two of H) for N ■ 14, in the worst case 
4xl4«56 independent equations are necessary, requiring N-m+1 points for the 
contour of each for applying the four successive conditions of the fields. 

In practice the symmetrical dlffusor is divided into N-m+-l touches (Figure 1). 

The solution of the systems of equations referred to allows the "forward 
scattering function" to be found (or advance diffusion) j. (o) along the 
two directions I and II shown in Figure 1. As already seen for the spherical 
drops, these functions are Indispensable for the calculation of attenuation. 

The values for j.,(o) calculated at F.U.B. at the frequencies presented 
in Table 3 and successive ones are presented in Tables 21-39 as a function of 
the various diameters of the spherical drops equivalent to the ellipsoidal ones 
with semiaxes a and b and eccentricity e . 

4. Specific Phase Displacement and Attenuation Along the Principal Planes 
The transition from individual scatter functions to attenuation and 
phase displacement, due to a medium composed of a discrete number of diffusors, 
is not immediate. As a matter of fact, the bonds between the fields and 
material defined by the equation 3 « c f , B • uH, t are immediate in 

the case of "continuous" media and not granular as in the present case. 
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Tlie difficulty has been overcome with the criterion of the equivalent 
of V. Dc Hulst [B.9], in which the equivalent is postulated between a 
dielectric layer and a layer of Che same thickness of spherical diffusors. 

Such a layer is presumed to be indefinite and transverse to the direction 
of propagation of a uniform plane wave [B.IO]. Tlie problem consists of 
determining the equivalent dielectric constant to be attributed to the di- 
electric, so that the mean radiation toward an observation point outside of 
Che layer (Figure 2) proves to be equal Co Chat which would exist at the same 
point with a real medium, Tlie solution of the equivalence is given by [B.12]: 

m*1 *JSfo)2»NKo 

where m is the complex refraction index of the medium (almost equal to 1). 

S(o) is the "forward scattering function" (presented in Tables 21-39 for the 
two directions of interest to us), is the propagation constant in free 
space. From this there is a direct transition to the propagation constant 
of a medium composed of a set of spherical diffusors given by: 

. I S.1CO) N (Di) 

K.0 i 

in the real case of ellipsoidal drops, assuming all of the drops to have 
rotational axes in Che transverse plane, defining the plane yz as the first 
principal plane (I) and the plane xy as the second principal plane (II) 

(Figure 1), Thus we can apply (S) along the two principal planes: (the 
index of the wave transmitted normally on the transverse plane, 9 = 90“ ) 



The application of (6) is based on an equivalence between the volumes of 
the ellipsoidal and spherical drops represented by means of [B.7]: 
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< ^ j! / 

- R a “ “S' ” 

IT ' * ’ 'V/ Cl* *'**'■ *'^ 

tJ ' <> • 9 

The values calculated for *H,i and ' j j ’ by 

means of j. presented In Tables 21-39 and the N(D^) of Table 2 are pre- 
sented in Tables 3-20. Ray's refraction index at the temperature of 20®C 
is also used in this case, 

5. Propagation of Radiowaves in Generic Polarization 

Tlie hypothesis of equioriented drops inclined with respect to the 
verticals, suggested by depolarization measurements, introduces a need to 
consider an inclination angle parameter " ^ " in the model propagating the 
electromagnetic waves [B.IO]: 

Since the state of polarization of an electromagnetic wave can be 
completely described by means of two special parameters a " and " <p ” 

[B.ll], the need arises to represent the propagating channel in a complete 
way, either from the physical viewpoint by means of the parame ter ‘t> and the 
propagation constants or from the electromagnetic viewpoint by means 

of the parameters "a " and "tj) 

The complete representation as a function of the parameters referred to 
is given by the matrix expression [B.ll]: 
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In which Yjj j- are the propagation constants described above, e is the 
connection length in Km, {|ui| is the unit matrix, lja{l is a matrix taking 
into account the state of polarization of the wave and of the angle of 

inclination of the ellipsoids, the terras of which a-e expressed by [B,ll]{ 

a,, * ' * 5 . 3 .. 0 . y i 1 -'33 ii» ‘ 

a, , 1 .9 r- ■ a ? tf Ei.. . ‘ . y 'ilr 

ft 

Table 40 represents the values to be attributed to 0 and <J» to obtain 
some special types of polarization of the electromagnetic waves which are of 
interest in practice, 

For the case of rectilinear polarizations, the values calculated by (7), 
with the hypothesis that ‘S> “ 0 , coincide with the values given in Tables 3-11 
for attenuation and Tables 12-20 for the phases. However, we can affirm 
that if the angle is limited to a few degrees, the values found for j. 
and ,, can now ba conaidarad valid. 

Considering circular polarization, the values of a and ii< in Table 40 
introduced in (7) tend to affirm that the attenuation and phase displacements 
in left circular polarization are equal to those in right circular polarization, 
and that both polarizations are independent of the angle of inclination of 
the drops. 

It is known from (7) that the Viilue of the attenuation and of the phase 
displacement in circular polarization can be conveniently expressed as the 
sum of a mean contribution (equal to the arithmetic mean of attenuation and 
phase displacements for 4>=o ), and a variation given by the model and 
argument of: 
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A * A 

The values of the mean A^j * II I‘ of attenuation and phase 

2 

displacement ate presented in Tables 3~20 at the frequencies previously 
given consideration. 

6. Propagation with Inclined Axes 

The case of earth-space propagation of electromagnetic waves in recti- 
linear or circular polarization can be treated in analogy with what has 
been said in the previous paragraphs . 

In the case of terrestrial radioelectrical connections, the direction 
of incidence of the transmitted wave is normal to the transverse plane 
(Figure 1; 0«90°) and, if we wish to consider inclined connections, it is 
sufficient to vary the angle of incidence 0 to obtain the values of the 
propagation constants y jj- ® » through the calculation program available 
in F.U.B.As an Indication Table 41 presents the values of -Yj) for 

three possible angles of inclination relative to waves transmitted in recti- 
linear polarization along the two principal planes of rain precipitation. 
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Table 1. Distribution of the Diameters of Drops as a Percentage of the 
Total Volume of Water (Laws - Parsons) 
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Table 7. Values of Attenuation (dB/Kra) at Precipitation Intensity 
R= 12.50(Tnm/h) 
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Table 8. Values of Attenuation (dB/Kra) at Precipitation Intensity 
R= 25.00(mm/h) 


21 



tP4^p..«p«CStC\iC^C^CSiCVfWfVCVCVCSiCN.Nrw 

ac»c>as»ooooc»<9«»ciao«»»;i»j» 










fV 

fS. 

»N 


fw 

rw 

Pw 


N 

fw 

N 

A. 


rs 


o 

o 

3 

3 

3 

3 

3 

9 

3 

3 

3 

3 

3 

3 

3 

9 

3 

3 

3 


« 

♦ 

• 

* 

« 

♦ 

• 

« 


♦ 

• 

# 

• 

• 

♦ 

♦ 

• 

♦ 

# 



▼ 



O 

> 


•0 

> 

-5 



3 

X 



A 


-> 


o 



X 

3 

3 



< 

% 

3 


(? 

X 

n 


CN» 

•S 


u 

1^ 


> 

W 

(V 


O 

**> 

tr 


> 

I\ 

O 

*7 

PS 

3 

?M 



-< 

▼ 


•3 

fv 

T' 




A 

C 

C 

PS 

PS 

V 

£ 

o 

> 

O 

> 


«-4 


A 

fN 

> 


•4 

ih4 

«*4 




**^ 

*•4 


««• 





30 a 009 '^» 9003000 -i» 9 S 9 
^ f' 4 ^^rv^T*-.«*vCV.C' 0 C»»O^fv,»wihi*^ 



»-« 


•• 

md 


rs. 

Pv 

cs 

PW 

fs 

-V 

<N. 

a; 

CS, 

fV P^ 

fV A. 

A. 


a 

a 

3 


3 

a 

s 

a 

a 

a 

a 

a 

3 

3 

a a 

3 3 

3 


♦ 

« 

» 

* 

# 

♦ 

♦ 

♦ 

# 

♦ 

♦ 

m 

♦ 

* 

« # 

• # 

* 


A 

*v 

> 

X 

*7 

0* 


o 

•o 

•o 


X* 

3 


> fo 

AO <o 



*s 



3 

-7 

3 


3 

?s» 

rw 

*A 

'T 

lA 

> 

PS 3 

lA X 

9 


■s. 

<c 

w 

■s 

3 

r> 

a 


A 

T 


rs 

M 

3 

3 V 

o o 

> 

< 

CN. 

X 

A 

▼ 

A 

3 



▼ 

iA 



p*. 

fs 

X X 

X X 

X 


^ <N4 

▼ 

O 

3 




-< 




•>4 






4 


^««^*ncv:vi:n*cs#^w NtCwrv^j^jesjiXi^vcuv'^ 
aoa 3*39 3 a'393'3 • 339 393 


o**w'>-^*>.o’rPwV* 5 AAA c 3 .«*!^riao 
— s:^»ircv< 9 «o<vi^*^«r^p^-^v 33 A 

^ A^>OGorv^i»*\<rvrvGOOOi 9 ^ 7 ^c%o ^3 


30339 <3 9 C 3300000^^^3 03 


oAoiroM^ oAotf\sir oaoaoao 
OOAAO 0 *^«^ 9«>»0 


Table 9. Values of Attenuation (dB/Km) at Precipitation Intensity 
R=? 50.00(ram/h) 
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Table 10. Values of Attenuation (dB/Km) at Precipitation Intensity 
R= 100 . 00 (mm/h) 
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at Precipitation Intensity 
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Intensity R= .25(mm/h) 
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Table 14. Values of Phase Displacement 


Intensity R= 2.50(mra/h) 


(® /Km) at the Precipitation 
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Table 15. Values of Phase Displacement ( “ /Km) at the I'recipitation 
Intensity R= 5.00(mm/h) 
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Values of Phase Displacement ( “/Km) at the Precipitation 
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Table 18. Values of Phase Displacement ( “/Km) at the Precipitation 
Intensity R= 50.00(mm/h) 
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Table 19. Values of Phase Displacement (“ /Km) at the Precipitation 
Intensity R= 100.00(mm/h) 
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Table 20. Values of Phase Displacement ( “ /Km) at the Precipitation 
Intensity R= 150.00(mm/h) 
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Table 21. Values of Forward Scattering Function at Frequency of 10 GHz. 
Key: 1-Diameter 
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Table 22. Values of Forward Scattering Function at Frequency of 15 GHz. 
Key: 1-Diaraeter 
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Table 24. Values of Forward Scattering Function at Frequency of 25 GHz. 
Key: 1-Diameter 
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Table 25. Values of Forward Scattering Function at Frequency of 30 GHz. 
Key: 1-Diameter 




38 



0 

1 



Table 26. Values of Forward Scattering Function at Frequency of 35 GHz. 
Key: 1-Diameter 
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Table 27. Values of Forward Scattering Function at Frequency of 40 GHz. 
Key: 1-Diameter 


40 





o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

9 

o 


o 

o 

o 

o 

o 

3 

o 

o 

o 

o 

o 


o 


• 

♦ 



• 


• 

« 

♦ 

♦ 

4 

♦ 

4 


W 

UJ 

w 

Ui 

ai 

Ul 


Xi 


XI 

XI 


yj 




o 


• 4 


» 

o 


*# 

94 




O 

o 

« 


> 


ir< 

A 

> 

*9 

'n 

9 

a 


r- 


14 

» 


A 

4^ 

4 

cc 







« 1 


4 


4 t 

4 

««» 

4 

!■« 

J 

o 




r4 


l>« 


f** 



o 

• 



«• 


« 

o 

« 

' 9 


9 

•1 

4 

4 


t 

« 



f* 

rv 

O 

0^ 






f** 

4 

f- 

«4 

V 

r% 


c 

tf. 

m 

«*• 

X 


% 

o 

4 

*•> 

Cl 


•n4 





<n 

.r 


4^ 

4> 


»• 

•O 


• 

• 

• 

• 

« 

• 

• 

• 

• 

• 

m 

• 

• 



o 

c 



o 

o 

o 

o 

o 

o 

o 

o 


") 


*» 


-3 


T> 

n 

~> 



-5 

”> 




o 


O 

« 


0-n 

mm 



-X 



o 


o 


o 

o 

o 

a 

a 

o 

o 

o 

o 


1 

1 

♦ 

* 

• 

« 

• 

4 

• 


9 

9 

9 


a 

111 

U.* 

u* 

it- 

lU 

u 

ir 

UJ 

•*l 

w 

u 

U 




ri<'« 


* • 

•* 

4 

c* 

> 

4\ 

•J 

l-» 

O 





•x> 

J 

o 

O 



’*4 

n 

r4 

m9 



tm ^ 

W 

O 



f*» 

•<i 1 

a 

.7 

4 

r4 

ry 


•c 

« 

;r 

vT- 

o 


f 


C 

O 


X* 

Ti 


r«* 



3 

w 

0 


■T' 

> 

•X 

4 

o 



9 


0^ 


iT 

f> 



y 

ir. 

h* 


4 




P f 

rk 

4 

•T 



**■ 

P-# 



rw 



* 1 

•TV 

O' 

3 

f t 

4J 



f* * 


■C 

* O 



■r 




- n 


^4 


♦ 

* 

4- 



O 

o 


o 

W' 





o 


o 




Table 28. Values of Forward Scattering Function at Frequency of 45 GHz. 


Key : 1-Diame ter 
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Table 29. Values of Forward Scattering Function at Frequency of 50 GHz. 


Key : 1-Diame ter 
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Table 30. Values of Forward Scattering Function at Frequency of 55 GHz. 
Key: 1-Diaraeter 
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Table 31. Values of Forward Scattering Function at Frequency of 60 GHz. 
Key: 1-Diaraeter 
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Table 32. Values of Forward Scattering Function at Frequency of 65 GHz. 
Key : 1-Diame ter 
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Table 33. Values of Forward Scattering Function at Frequency of 70 GHz. 
Key: 1-Diameter 
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Table 34. Values of Forward Scattering Function at Frequency of 75 GHz. 
Key; 1-Diameter 
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Table 35. Values o£ forward Scabterlng Function at Ftaguency o£ 80 GHz 
Key : 1-Diame ter 
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Table 36. Values of Forward Scattering Function at Frequency of 85 GHz 
Key: 1-Diameter 
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Table 37. Values of Forward Scattering Function at Frequency of 90 GHz. 
Key: 1-Diaraeter 
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Table 38. Values of Forward Scattering Function at Frequency of 95 GHz. 

V 

Key; 1-Diameter 
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Table 39. Values of Forward Scattering Function at Frequency of 100 GHz. 
Key: 1-Diaraeter 
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Table 40. Values of the Parameters tp and a for Various Types of Polarization 
Key: 1-Linear horizontal, 2-Linear, bisector I and III quadrant, 3-Circular 

left, 4-Linear, bisector II and IV quadrant, 5-Circular right, 6-Linear vertical. 
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